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Abstract 

We consider the 2D SU{N) principal chiral model and discuss a vortex condensation 
mechanism which could explain the existence of a non-zero mass gap at arbitrarily small 
values of the coupling constant. The mechanism is an analogue of the vortex condensation 
mechanism of confinement in 4D non-Abelian gauge theories. We formulate a sufficient 
condition for the mass gap to be non- vanishing in terms of the behaviour of the vortex free 
energy. The SU(2) model is studied in detail. In one dimension we calculate the vortex 
free energy exactly. An effective model for the center variables of the spin configurations of 
the 2D SU(2) model is proposed and the Z{2) correlation function is derived in this model. 
We define a Z(2) mass gap in both the full and effective model and argue that they should 
coincide whenever the genuine mass gap is non-zero. We show via Monte-Carlo simulations of 
the SU{2) model that the Z(2) mass gap reproduces the full mass gap with perfect accuracy. 
We also test this mechanism in the positive link model which is an analogue of the positive 
plaquette model in gauge theories and find excellent agreement between the full and the Z(2) 
mass gap. 
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1 Introduction: Motivation and Problems 



The nature of a non-zero string tension (ST) at any value of the bare coupling constant is one of 
the most striking puzzles of QCD. An analogous question exists in two-dimensional (2D) non- 
Abelian spin models with continuous symmetry group where the mass gap (MG) is expected to 
be non-zero at any temperature. Despite huge efforts in both solution of this problem 

has not yet been found. A promising conjecture which could explain the existence of a non- 
zero ST in QCD and thus permanent confinement of quarks is based on a vortex condensation 
mechanism. In this article we adjust the lattice formulation of this mechanism Q to the case 
of 2D non-Abelian spin models. Though this paper is in a spirit of [g], our treatment of the 
problem differs in some aspects. We propose a simple effective model for the Z{2) degrees of 
freedom in the SU(2) principal chiral model which explains how this mechanism could work in 
the weak coupling region of the lattice model. 

The idea that condensation of vortices may be responsible for confinement of static quarks in 
non-Abelian gauge theories with non-trivial center appeared already in the late seventies |3], Q. 
The essential concept was taken from Z(N) lattice gauge theory (LGT) where the Wilson loop 
was known to obey an area law in the strong coupling region. Since Z(N) forms the center of the 
SU(N) group it was suggested that Z(N) vortices can be also present in the more complicated 
SU(N) theory and play an essential role in generating a non-zero ST. However, Z(N) LGT 
undergoes a phase transition at weak coupling to the deconfinement phase with the Wilson loop 
obeying a perimeter law. Since Z(N) vortices in Z(N) LGT may have a thickness of only one 
lattice spacing its contribution to the free energy becomes negligible in the weak coupling region 
where the system is well ordered at small distances. What concerns SU(N) LGT, attempts to 
calculate an effective Z(N) theory at small coupling making a perturbative expansion around 
Z(N) solutions of the Yang- Mills equations only lowered the critical coupling but did not remove 
it to zero as is expected to be the case for the correct confinement mechanism Q. 

In theories with continuous symmetry group like SU(N), vortices may have, however, a 
thickness of not only one but many lattice spacings. Hence, there is the possibility to generalize 
the naive mechanism of Z(N) models and include all the possible vortex configurations present 
in SU(N) theories. Such a theory of confinement was developed in [|| where also a theorem was 
proved which makes a link between the behaviour of the Wilson loop and the vortex condensate. 
In this theory the definition of a vortex is actually not important. The only important issue 
is a change of vorticity, and the vortex condensate is defined as the free energy of such a 
change introduced by special singular Z(N) gauge transformations. Over large distances typical 
configurations look like Z(N) vortices, i.e. the basic field variables separated by such a distance 
are rotated by a Z(N) element relatively to each other. The vortex condensation mechanism 
is manifestly gauge invariant and presumably gives a nice explanation of the coexistence of 
confinement at large scales and perturbative behaviour at short range. The mechanism was 
investigated in many papers and some results supporting its validity were found f5[- fl2]| : 

• A strong coupling expansion of the vortex free energy up to the 12-th order demonstrates 
that vortex configurations produce a ST which coincides with the full ST up to this order 
||. This result is gauge independent. 

• The Z(N) Wilson loop was shown to carry all the ST to all orders of the strong coupling 
expansion ||, at least in the electric gauge. It was argued that Z(N) Wilson loops in 
different gauges differ only by perimeter contributions. 
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If the magnetic-flux free energy vanishes in the limit of a large uniform dilatation of a 
torus, the vortex free energy always decreases exponentially. It is sufficient to produce 
confinement. Using this property one can rigorously prove the lower confining bound in 
three-dimensional U(N) LGT 0. 

Simple intuitive ideas as well as some analytical results show that a vortex mechanism is 
likely to lead to confinement at weak couplings also in the positive plaquette model and 
the Mack-Petkova model, which eliminates certain Z{N) magnetic monopoles in SU{N) 
LGT §. 

The technique to evaluate the contribution of vortices of arbitrary thickness to the expec- 
tation value of any observable was developed in H. Such a "preaveraged" Wilson loop, 
i.e. calculated solely on the vortex contributions, exhibits confining behaviour. While not 
rigorous, this result potentially refers also to the weak coupling limit of the 3-D SU{2) 
model. 

It has been proved that even the classical, though not naive limit of SU(N) LGT includes 
bare vortices in the continuum Lagrangian. They are labeled by the nontrivial center 
elements of SU(N) and are supported on closed 2D surfaces in four dimensions [|l0|]. A 
one-loop expansion in a particular background of such vortices shows instability of the 
vacuum implying that vortices must condense (become "fatter") in the quantum theory 
already at two loop order. 

Recently, the vortex condensation mechanism of confinement was studied in the so-called 



maximal center gauge of SU(2) LGT [11|. The authors of these papers claimed that the 
Wilson loop computed with the center projected Z(2) gauge field degrees of freedom carries 
almost the whole asymptotic ST. Consequently, excluding all Z(2) vortices identified after 
projection leads to a vanishing ST. 

The SU(2) partition function can be rewritten in the form of coupled SU(2)/Z(2) and 
Z(2) models which allows to give a proper interpretation of different Z{2) excitations in 



the original model [12|. Using plausible assumptions one can establish a link between these 
excitations and the behaviour of the sign of the trace of the Wilson loop. It was shown 
via Monte-Carlo (MC) simulations that the sign of the trace of the Wilson loop carries all 
the information about the asymptotic behaviour of the fundamental string tension. 

It should be stressed that earlier papers on the vortex condensation mechanism of confine- 
ment in continuum QCD cannot be regarded as an explanation of confinement ||13[| : In absence 
of a non-perturbative definition the introduction of vortex configurations into the QCD La- 
grangian seems to be a completely ad-hoc procedure. Thus such models are not able to give an 
explanation in terms of dynamical reasons why vortices should become fat, i.e. why they are 
condensed. 

In general, there are two ways of looking at the vortex condensation mechanism in lattice 
QCD. The first one results from the desire to find an analogy with the continuum theory and 
interprets the vortices responsible for confinement as an analogue of the Nielsen-Olesen vortices 
and the QCD vacuum as the so-called spaghetti vacuum. The second approach is based on the 
similarities between lattice QCD and 2D non-Abelian spin models where, while not so close to 
the continuum, one can give precise mathematical definitions to all quantities involved. Following 



3 



the latter idea, Mack and Petkova |l|] formulated a condition which could be called confinement 
mechanism by a vortex condensate. We would like to emphasize that while the relation of this 
approach to the spaghetti vacuum is rather vague at the present stage of affairs the connection 
to 2D spin models is straightforward since one has precise definitions in both cases. Moreover, 
a quantity like the vortex condensate is expected to be a genuine non-perturbative quantity 
and thus it should be clear that it can be given a precise meaning only in a non-perturbative 
approach such as the formulation of QCD on a lattice. 

In this article we consider the vortex mechanism in some details on the example of 2D spin 
models. The reason to deal with these models is the following: The nature of the MG is unknown 
despite the claim in the literature that its exact physical value is known. A vortex condensation 
mechanism is one possible candidate to explain the phenomenon of a non-zero MG. Moreover, 
this mechanism has an analogue in gauge theories. In view of all similarities between 2D spin 
models and AD gauge theories we think it is useful and instructive to study this mechanism on 
the example of simpler 2D models. 

This paper is organized as follows. As a first point we give a definition of the vortex free 
energy in terms of the original spin configurations. We introduce a vortex container which 
has the topology of a ring in 2D and which is specified by certain boundary conditions (BC). 
We prove a sufficient condition for the correlation function to decrease exponentially which is 
precisely the analogue of the Mack-Petkova theorem in LGT. This is done is Section 2. 

Having identified the exponential decay of the vortex free energy as a sufficient condition for 
producing a non- vanishing MG the question arises what are the configurations of the spin field 



responsible for this rapid change of vorticity. Following the ideas of || 14] one can introduce a 
Z(N) correlation function which measures the effect of vortices on the full correlation function. 
In spin models with global symmetry this is conceptually easier since we do not have to fix 
a gauge to uniquely define such a quantity. One can give the corresponding arguments [14 
that this Z(N) correlation function defines completely the large distance behaviour of the full 
correlation if a vortex mechanism is responsible for the non- vanishing MG. We shall construct 
such a correlation and derive an effective Ising-like model for Z(2) excitations in Section 3. 
Furthermore, we calculate the correlation function numerically for the case of the SU (2) principal 
chiral model using MC simulations. We find that the MG extracted from the Z{2) correlation 
function agrees almost completely with the asymptotic MG of the 577(2) model. A more delicate 
question is what one can expect in the case of the positive link (PL) model where all the thin 
vortices are eliminated. We introduce this model in analogy to the PP model || and study it 
in the SU{2) case. In particular, we define and calculate Z{2) correlations in the same way as 
in the full 577(2) model. All these questions are subject of Section 4. 
Our conclusions and discussion are presented in Section 5. 



2 Vortex mechanism in spin models 

In this section we give a precise formulation of the vortex condensation mechanism in two- 
dimensional spin models. A first hint that such a mechanism could be crucial for the existence 



of a non-zero MG came from the famous paper by Dobrushin and Shlosman [15]. These authors 
pointed out that the Mermin- Wagner theorem on the absence of spontaneous magnetization in 
two-dimensional spin systems follows from the intuitive idea that in such systems long Peierls 
contours cost only little free energy by making them thick. This is possible because the spins 
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can rotate slowly due to the continuous nature of the symmetry group. 

In the following we derive a sufficient condition for the MG in two-dimensional SU (N) 
spin models to be non-vanishing. Our derivation follows closely the one of the corresponding 
condition in LGT by Mack and Petkova 0. Nevertheless, we choose to adduce it here because 
of pedagogical reasons and to make this paper self-contained. 

For definiteness we consider the SU(N) principal chiral model on a two-dimensional periodic 
lattice A D . Its partition function is given by 

Z((3)= [ [J A"(^x0exp[/?^ReTr(C4.[4 +n )], U x G SU(N). (1) 

J x£A D I 

/3 is the coupling constant, Dfi denotes the normalized Haar measure on SU(N) and the inte- 
gration extends over all sites x of the lattice. The sum in the exponent is a sum over all links 
/ = (x,n) of A D . According to the conventional scenario, the fundamental correlation function 
of spins separated by a distance R 

r K (/?) = (ReTr(C7 4)) ( 2 ) 

is expected to decrease exponentially at any value of the coupling constant /3 if R is sufficiently 
large 

T R (p) ~ exp(-m c (/5)i?), (3) 

with m c ((3) the MG. 

Let us consider now sets of links on A D which form closed loops on the dual lattice as shown 
m fig. g. Following we call the region of the original lattice enclosed by two such loops a 
vortex container T which has the topology of a ring in two dimensions. We arrange many such 
containers between the lattice sites and R. Different containers may touch but not intersect 
each other. Let dTi be the boundary (the set of circles • and squares ■ in fig. [l]) and Aj the 
interior (the set of crosses + in fig. |l|) of the i th container Tj. We define the complement A c 
of A^* as A c = A- D /]^[ i Aj. To make the following derivation more transparent we rename the 
SU(N) variables U x . Spins belonging to one of the vortex containers Tj are renamed U' x , spins 
lying in the complement A c U x . 

We start by rewriting the path integral expression for the fundamental correlation function 
O) in terms of the variables UL and U r 



Tr(P) = \ f II Dfi{U x )ReTT(U Ul) exp 



f3 ]T ReTr(U x ul 



ieA c 



n / n ^e'iwwmu ii *mo • w 



The group 5-function is necessary to avoid double integration over spins defined on the boundary 
dTi of one of the containers. The product runs over all containers arranged between the lattice 
sites and R. In the sum ^ in the exponent of the inner integral, links in the boundary of the 
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Figure 1: An example of a vortex container T on a two-dimensional lattice. The vortex container 
is the set of all lattice sites which are enclosed by the two closed loops depicted as dashed lines, 
i.e. the set of all crosses +, circles • and squares ■. The interior A of the container T is the set 
of all crosses +, whereas the boundary dT consists of all circles • (= the outer boundary dT° ut ) 
and squares ■ (= the inner boundary dT m ). 



container Tj are omitted. The inner integral equals the partition function Z(Ti, U x ) defined on 
the container Tj with boundary conditions U x , x G dT{, 



Z(T i ,U x )= I H Dii{U' x ) exp 



/3^'ReTr(£/^ +n ) 



J] S^U- 1 ). (5) 

xedTi 



An important property of Z(Ti, U x ) is the invariance under SU(N) transformations of the bound- 
ary conditions. In particular, Z(Ti,U x ) remains unchanged under the transformation 

U x -> uj- 1 U x , x G dTi, (6) 

with u>i an element of the center Z(N) of SU(N). Let us perform now the following variable 
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substitution 



U, 



N c 

k=i 
N c 

U x JJ UJ k 

k=l 



if 



x < Ai 



Aat < X. 



(7) 



The notation Aj_i < x < Aj means that the lattice site x belongs to the region which is enclosed 
by the interior Aj_i of the container T_i and the interior Aj of the container Tj. N c is the 
total number of containers. In the special case of U$ and Or, the substitution ([?]) results in 
Uq = Uq Yli u)i and Ur = Ur. For the correlation function (|j) we can write 



Tr{P) = \ / II D^U x )ReTr(U U R )exp[P^2 ReTr ( U M + n)} H^Z^U: 



(8) 



xeA c /gA c 
with boundary conditions U x on the container Tj 

N c 

V% = U x 11 cu k if 

k=i+l 

N c 
k=i+l 



G dT" 



G dT 



■out 



(9) 



0T? n and dT° ut are the inner (the set of squares ■ in fig. [l]) and outer boundary (the set of circles 
• in fig. ||) of the container T respectively. We further simplify the boundary conditions (||) by 
applying the transformation U x — > U x Ui Y\ k=i+1 oj^ 1 under which Z(Ti, U£) remains unchanged. 
Thus we finally can write 



U x LJi 

u x 



if 



x e dT? 
x g dT 



out 



(10) 



for the boundary conditions XJ£ of the container Tj. Since Wj G Z(N) is arbitrary, we may 
sum (integrate) over Wj in (h) using the normalized Haar measure on Z(N). Using the trivial 
relations 



ReTr(U U R ) < Tr(l 



< max 



E^z(t u u-) 



(11) 



and with the help of the identity 

Z / ]jD^U x )eM0j2 Re ^ U - U Un)}UJ2 Z ^ U - 



(12) 



ieA c 
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we derive the following bound for the correlation function (Q) 



11*00)1 < Tr(l)[]max 



(13) 



The maximum on the right hand side of (O) is to be understood as the maximum under the 



boundary conditions (1C). In order to make this result more transparent we will specify it for 



the case of SU(2). Then u>i can take values ±1 and the bound (13) can be written 

\r R ((3)\ < Tr(l) 17 max |^(T i5 C/^)|, (14) 

Ux 

where we introduced 

Vi(Ti, U") = . (15) 

-L "T q i 

qi is called the vortex free energy and is given by the ratio of partition functions defined on the 
vortex container Tj with boundary conditions U£ ([l0|) 



Z(T h U: 



(16) 



The physical meaning of % (and more generally of Vi) is rather obvious: qi measures the change 
of free energy exp(— AF) introduced by a singular transformation U x — » U x Ui on one of the 
boundaries of the container Tj. Let Li be the diameter of the container Tj. Then we say that 
the transformation U x — ► U x Ui performed on one of the boundaries of the container introduces 
a vortex of thickness Lj into the system and qi measures the energy needed to create such a 
vortex. 

We are now ready to explain what we call vortex condensation mechanism for generating a 
non-zero MG in 2D SU(N) spin models. Let l/ max be the maximum of \Vi\ under the boundary 
conditions U x defined in (|lOD; i.e. V™"^ = max^ \ Vi\. Suppose that there exists such a mass 
m v that for each container Tj and for sufficiently large diameter Li, the maximum V^ nax behaves 
according to 

V™™ ~ e W (-m v ((3)Li). (17) 

Then, with R = ^ L«, the bound (|i~4] ) for the correlation function reads 

\T R (p)\ < const -exp(-m^ ((3)R), (18) 

which is the expected exponential decay. We term m v vortex MG to distinguish it from the 
genuine MG m c extracted from the correlation function. If m v = m c , condensation of vortices 
can be made responsible for generating a non-zero MG in 2D SU(N) spin models. 

The crucial quantity in a vortex condensation mechanism is the vortex free energy qi which 
is originally defined on the container Tj but which may be calculated on any lattice having the 
same topology as the container; i.e. the topology of a ring. Thus, we have to consider a lattice 
with periodic boundary conditions in one direction and fixed boundary conditions in the other 
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Figure 2: A lattice (periodically closed) which has the same topology as the container shown 
in fig. ||. Spins defined on circles • and squares ■ are fixed to given values. In the case of the 
vortex condensation mechanism introduced in ref. 0, the coupling (3 is changed to —j3 on the 
depicted links. 



direction. A corresponding lattice is shown in fig. |2[ The vortex free energy q on this lattice is 
given by 



Z(U x u>) 
Z(U X ) 



(19) 



where in the denominator values U x are assigned to spins depicted as circles • and squares ■ in 
fig. ||, while in the numerator spins living either on circles or on squares are fixed to values U x uj, 
with uj 6 Z(2). The partition function Z is defined in (||). 

Finally, we think that still some comments are needed at this place: 

• Configurations which rotate slowly from one center element to another one on some char- 
acteristic length scale are termed "thick vortices". According to the paper by Dobrushin 



and Shlosman [15] they produce disordering effects which are sufficient to enforce the long 
distance correlation function to fall off to zero at any coupling; i.e. they guarantee the 
absence of magnetization. However, a priori there is no reason to believe that thick vor- 
tices provide an exponential fall off. It might be that for (3 — > oo these configurations can 
only account for a power law decay of the correlation function. Two conclusions can be 
drawn from such a scenario. 1) Thick vortices are not responsible for a non-zero MG at 
arbitrarily large values of (3 or 2) they are at small (3 but the system undergoes a phase 
transition to a massless phase at a finite value of (3. 

The above vortex condensation theory does not give an exact definition of a vortex. The 
only important issue is a change of vorticity over some characteristic length scale. To 
specify a vortex completely, one has to define the rate of this change and, possibly, the 
dependence of the characteristic length scale on the bare coupling (3. It is reasonable 
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to assume that on this length scale, the SU(N) spin model looks like an effective Z{N) 
model for special Z(N) excitations which one should be able to extract from the original 
configurations. 

2.1 ID model 

Before finishing this section we will demonstrate on a simple example how the mechanism de- 
scribed above works in practice. The simplest example we can imagine is the SU(2) spin model 
in one dimension where it is well known that the correlation function shows an exponential fall 
off at any value of the coupling constant (3 and the corresponding MG is given by 

with I n modified Bessel functions. In one dimension a vortex container is simply a chain of spins 
with Dirichlet-like boundary conditions. More precisely: let L be the length of the ID chain 
and let us fix the spins on the boundary to some arbitrary value 

U x=0 = W x , U X=L = W 2 . (21) 
Define W = WxW%- On a finite lattice, the partition function Z can be calculated exactly 

oo 

Z(W) = P~ L (.2n + l)xn(W)[I 2n+1 (P)} L , (22) 

71=0,1/2,... 

where XniW) ls the character of the n th representation of SU{2). Perform now a nontrivial 
Z(2) transformation on one of the boundary spins of our chain, e.g. W 2 — > ojW 2 = —W 2 . The 
character XniW) transforms according to 

Xn(W) -> Xn(-W) = (-lf n X n(W). (23) 

For the vortex free energy q ( jlBf ) we obtain 

9 ^ K1 _ 4X1/2W (||)%..., 

where the terms which vanish faster in the thermodynamic limit have been neglected. The 
maximum of V fll4| ) occurs at W = 1 and we finally find 

m&xV ~ exp[— Lin 
w 

Thus, in the ID model the change of V introduced by the creation of a thick vortex shows an 
exponential decrease with a MG which equals the MG extracted from the correlation function 
(|20|). Of course, in one dimension this result is in a sense trivial but it shows nevertheless 
transparently how the idea described above works. 



h(P) 
W) 



(25) 
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3 Effective model for Z{2) excitations 



In this section we address the question which are the spin field configurations playing a crucial 
role in a vortex condensation theory. In the case of gauge theories, it was suggested long ago 
that certain Z{N) excitations of the gauge field are responsible for an area law behaviour of the 
Wilson loop [14}| . Here, we present an effective model for Z{2) excitations in the 2D SU(2) spin 
model and establish a link to the vortex condensation mechanism. The main assumption of this 
model is supported by MC data which will be presented in next Section. 

As a first step we fix proper boundary conditions of the two-dimensional lattice. It is con- 
venient to take a periodic lattice in y-direction and to fix the SU{2) spins to elements of the 
Z{2) subgroup if x = and x = L. L is the linear extension of the lattice. It should be stressed 
that fixing boundary conditions is not necessary for systems with global symmetry since the 
Mermin- Wagner theorem guarantees independence of the results on BC in the infinite volume 
limit. In general, a dependence on BC should vanish faster than any long-distance correlation 
function. In our case fixing boundary conditions is only a question of convenience and proper 
definitions. 

We start by rewriting the partition function of the SU(N = 2) principal chiral model ([l]) 
using the representation 

U x = z x U x (26) 

for the spin field U x £ SU(2), where z x is a Z(2) element and U x £ SU{2)/Z{2). In the new 
variables the BC for z x are free ones while for U x we have to impose Dirichlet BC in x-direction. 
For the invariant measure on the SU{2) group we write 

DKU X ) = \ D »(U X ), (27) 

{z x }=±l 

where D[i(U x ) is an invariant measure on the SO(3) group. Let us recall that the invariant 
measure on the SU{N)/Z{N) group coincides with the SU(N) measure up to the restriction 

~ < arg(Triy < ^. (28) 

In other words, in the invariant [/^-integration the trace of the fundamental characters is re- 
stricted to positive values. Performing now the summation over the Z(2) elements z x one can 
rewrite the partition function (|l|) as (up to an irrelevant constant) 

Z W) = E E f\{D^U x ) e ^zZ s ^^+n)\X\ s i- ( 29 ) 

C is a set of closed loops and we have introduced a new Z(2) link variable si. This representation 
for Z{(3) (and a similar one in the case of SU(2) gauge theory) was derived and investigated 
in |l6| . In terms of the link variable si the singular transformations discussed in section [2] may 
be defined as a change of the sign of si on the boundaries of a vortex container T (only in 
the action). The sum over C in (p9|) is a sum over all closed loops (including all their possible 
products) taken with an appropriate weight. It is defined exactly in the same way as in the two- 
dimensional Ising model where much is known about the properties of such a loop expansion. 
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However, unlike the Ising model, the coefficients defined on the links of closed loops are not 
constant. Moreover, if we consider the model (29) as an Ising-like model with a fluctuating 
coupling constant, we observe that the coupling is not positive definite (despite TiU x > this 
is not the case for the character of the product of two group elements). At least at first glance, 
this could imply that fluctuations of the link variable S[ may persist down to weak coupling and 
cause the disorder which is needed for an exponential behaviour of the correlation function. In 
fact, this is precisely what we are going to work out. 

As a first step we perform the sum over the link variables sj. It results in 

z(p) =e / n^/ i (^)n c ° sh ^ Tr (^^+n)]n tanh [/ 3 ^(^^V)]- ( 3o ) 

c J x i iec 

For sufficiently large values of (3 we obtain 

Z{(3) = ^(tanh2/?)l £ l I \[D l i(p x )J[ew{P\^(pxUl^)\]'Y[<ri + 0(e' 4/3 ), (31) 
c •> x i iec 

where 

^ = sign[TY(^C?] + J]. (32) 

iFiom the last equations we see that at large values of (3 the original partition function, and so 
the free energy, can be written as a product of two partition functions - the partition function 
of a SU(2)/Z(2) model and the partition function of an Ising-like model. Precisely 

Z(J3) = z^/^Z 1 , (33) 

with 

Z sumim = f ]jD»(U x )]Jex V [P\nU*ul + n)\] (34) 

J x l 

and the Ising-like partition function 

Z 1 = ^(tanh2/?)l £ lF(£). (35) 
c 

F(C) is defined as an expectation value 

F(£) = (]Jai}su(2)/z(2) (36) 
lec 



which has to be evaluated in the ensemble (|34|). In the case of the fundamental correlation 
function (0) for iV = 2 we follow the same strategy and obtain for large (3- values 



T R (P) = (Z su ^ z ^Z I r 1 Y,( t ^ h2 ^ Vl £( tanh2 ^) |£| fl[D»(U x )Tr(UoU 



i) >< 



eMPlMUxUl^nH^H^ + 0(e~^). (37) 
lev i&C 
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Here, the sum over V is a sum over all paths (including all their possible products) connecting the 
lattice sites and R. In addition, one has to sum up over all possible closed loops on the lattice 
which must not have a link in common with a given path V . Again, this is in full accordance 
with the expansion of the correlation function in the Ising model. To see this explicitly, let us 
write down the corresponding formulae for the Ising model. For the partition function one gets 
an expansion in terms of loops (up to a constant) 

Z Isin9 = ^(tanh/3)l £ l, (38) 
c 

while the correlation function reads 

T R sin9 (f3) = (Z 7 ^)-i^( tanh/3 )|P| (tanh/3)l £ l. (39) 

p c/iev 



These formulae have to be compared with the asymptotic expansions ( pl| ) and ([37]) in the case 
of the SU(2) model. In the Ising model it is known that for /3-values below (3 C the correlation 
function is well approximated by the expression 

r^(/3</3 c )«^(tanh/3)l^l, (40) 

Vc 

where we have to take into account only the so-called connected graphs V c . In the following we 
shall show that one can arrive at a similar expression for the correlation function in the SU (2) 
model by making suitable assumptions. Firstly, we introduce a Z{2) correlation function in the 
SU (2) ensemble according to 

T Z /\(3) = (z z R ) (41) 
with z x = sign(Tr£/a;) and claim that it reproduces the correct long distance behaviour of the 



full SU{2) model, i.e. the MG extracted from @ coincides with the full SU(2) MG. This 
assumption will be confirmed in the next section by numerical results. It leads to the following 
expression for the asymptotic expansion of the SU(2) correlation function (RTj) 



= ^ 7 E( tanh2 / 3 )' P| E (tanh2/3)l £ l ([[ a t ]J ai) su{2)/z{2) , (42) 
V c/iev lev lec 

where the expectation value is defined in the ensemble (Bl). Our next assumption is that for 



(3- values smaller than j3 c (42) can be written as a sum over connected graphs like in the Ising 
model 

Fr(P < p c ) « ^(tanh 2/3)1^1 ( ]J a t ) su(2)/z{2) . (43) 
V c leVc 



This formula has to be compared with the corresponding formula (40) in the case of the Ising 
model. In the SU{2) case the conventional scenario means (3 C = oo. Thus we expect that 
approximation (^) works rather well in the whole region of the bare coupling. Of course, its 
goodness is determined by the behaviour of the expectation value on the right hand side of (^) . 
Suppose for a while that (HzeP a i)su(2)/z(2) is small enough to cancel the fast growing number 
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of connected paths with length \V C \ and to guarantee a fast convergence on the right-hand side 
of (f43[). Then, it is easy to give a prediction for the MG. To leading order one finds for j3 < j3 c 



(tanh2/?) fl ( JJ o\ )su(2)/z(2) 

in 



(44) 



where V m i n is the shortest path between the sites and R. Next to leading order corrections 
can be obtained from (^3|). 

Let us discuss now some of the issues involved with the above Ising-like effective model as 
well as some possible physical scenarios. 

1. A link of the above effective model to the vortex condensation theory presented in the 
previous section can be established as follows: In analogy to the Z(2) correlation function 
let us define the Z{2) vortex free energy. To do this, we rewrite the ratio of partition 
functions q defined in ( |i~9| ) as an expectation value of an appropriate operator which can 
be expanded in a sum over SU(2) representations. Making use of the decomposition (EG) 



(28) and proceeding along the same line as in the case of the correlation function we 



are able to express the vortex free energy in terms of expectation values defined in (pq) 



and (43). Thus, if it turns out that condensation of vortices is responsible for a non-zero 
MG and this MG can be extracted from (|37|) and (|43|) respectively, then the Z(2) degrees 
of freedom reproduce the exponential decay of the vortex free energy as well. However, 
this is only possible in a SU(2)/Z(2) background, since the Z(2) degrees of freedom alone 
cannot account for an exponential fall off. 

2. If Z(2) degrees of freedom indeed play a crucial role in generating a non-zero MG, this 
may have some strong impact to SU{2)/Z{2) models. Consider for example the following 
variant of an 50(3) model which is an analogue of the lattice gauge model introduced in 
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Z SO{3) = [U D ^ e ^/Z S ^ U ^+n)] = 

{s t }=±l J x I 

= [j[D»(U x )]Jco S h[(3Tr(U x U x+n )}. (45) 

J x I 

At large (5 we obtain 

z so(3) = f Y[D»(U X ) HeMP\Tr(uM + n)\]- (46) 

J x I 

Speculations about the coincidence of SO (3) and 577(2) models are based on the widely 
accepted belief that the continuum limit has to be taken at (5 — > oo, where the naive 
continuum limits of both models coincide. But since the 50(3) model ( [46|) lacks of Z{2) 
degrees of freedom, a possible non-zero MG cannot be explained with the above described 
effective model. Thus, either SU{2) and 50(3) models have different continuum limits, 
or the Ising-like part (^) of the SU(2) partition function ( |33| ) as well as the correlation 
function (f43|) must become trivial in this limit. The latter scenario means that the MG 
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extracted from the Z{2) degrees of freedom vanishes above some pseudo-critical value of 
(3 while the full SU(2) MG does not. This seems to be very unlikely, since at least for 
finite values of (3 the Z(2) correlation function carries the whole MG. This will be shown 
by numerical results in the next section. 

3. Suppose for a moment that the Z{2) correlation function indeed reproduces the full MG. 
Comparing formulae (^) and ( j38l) one can determine an effective Ising coupling in the 
SU (2) model. Two scenarios are possible. If the effective coupling is less than the critical 
coupling of the Ising model for arbitrarily large values of the bare coupling (3, the SU(2) 
model is always in a phase with a non-zero MG and the continuum limit may be taken 
at (3 — > oo according to the conventional scenario (at the same time the effective coupling 
has to approach the critical coupling of the Ising model, otherwise the very existence of 
a nontrivial continuum limit becomes problematic). It may happen, however, that at 
some large value of (3 the effective Ising coupling becomes larger than the critical coupling 
of the Ising model. Then above this /3-value the system is in a massless phase and the 
conventional scenario is broken since one has to realize the continuum limit at this finite 
/3-value corresponding to the critical effective Ising coupling. It is interesting to mention 
that in this case if it is possible to construct the massless continuum limit by driving the 
bare coupling to its critical value from above, one may expect that this limit coincides 
with the continuum limit of the SO(3) model. 



4 Monte-Carlo study of Z(2) and SU(2) mass gap 

In this section we will show by results of MC simulations that the long distance behaviour of 
the full SU{2) correlation function (^) coincides with the long distance behaviour of the Z{2) 
correlation function (|4l|). This is done by calculating the full SU(2) correlation length £,su(2) 
(which is the inverse of the MG m) and by comparing it with the correlation length £z(2) 
extracted from Z{2) degrees of freedom. As will be seen, both quantities - Csu(2) an d £,z(2) - 
agree within errorbars. This result confirms the assumption made in the previous section that 
the Z(2) degrees of freedom carry the full SU{2) MG. Furthermore, we consider the PL model 
and demonstrate that Z(2) excitations reproduce the full SU(2) MG in this model as well. 



4.1 Standard model 

To calculate the full SU (2) correlation length Csu(2) we follow the method presented in 
We parameterize our field variables U x £ SU(2) according to 

U x = u (x) + iu(x)a, u x = (u (x),u(x)), (47) 

^To be precise, in ref. the correlation length £ is defined as the second-moment correlation length and not 
as the inverse of the mass gap m which can be determined by fitting the falloff of the zero-momentum correlation 
function to the cosh-behaviour appropriate for time-periodicity. However, both quantities are expected to show 
the same scaling behaviour. Surprisingly, it is found empirically that the two definitions of £ do not only scale in 
the same way but agree within less than 1% |cj. In this article we employ the definition of £ presented in 
[|19j since it is less CPU-time consuming and we feel free to call it the inverse of the MG. 
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with a the Pauli matrices. We then determine the susceptibility XSU(2) according to 



XSU(2) = ^<G>>-) 2 > ( 4§ ) 



and the analogous quantity at the smallest non-zero momentum 
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F S U(2) - 



2nix 2 /L |2 



X 



}. (49) 



In terms of these two quantities the second-moment correlation length £su(2) is given by [19] 

ten = ( W^w-' V^. (50) 

\ 4 sin 7T/L / 

L is the linear size of the lattice and V = L 2 the number of sites in the lattice. In addition, we 
calculate the internal energy 

E S U(2) = T X (P) = (Tr(U ul)) (51) 

corresponding to the correlation function of spins separated by one lattice spacing. In the case 
of the Z(2) correlation length £z(2) we consider the same ensemble of SU(2) degrees of freedom 
U x but instead of calculating the quantities \ (0) an d F ( [49] ) with the parameters u x ( |47| ) we 
use the Z{2) degrees of freedom 

z x = siga(TrU x ) = sign(u (x)). (52) 



Cz(2) is then determined by inserting Xz{2) an d Fz(2) (pQI) and Ez(z) ^ s given by (zqZi). 

In order to simulate the 2D SU(2) principal chiral model defined in (|l|) with N = 2 we 
use Wolff's cluster algorithm |21|| . The model is considered at different values of the inverse 
coupling constant (3. We choose periodic boundary conditions. In all cases we start with a 
random configuration and apply at least 10 4 warm up sweeps. We generate 5- 10 5 configurations 
and measure the quantities of interest in every configuration since the cluster algorithm is known 
to show small autocorrelation times. To estimate errorbars we use the jackknife sub-ensemble 
analysis. In table Q we show the results for £577(2) an d £z(2) f° r several (3- values and lattice sizes 
L = 128 and L = 256. The data in the column denoted with is taken from and serves as 
reference results for Z,su(2)- As is clearly seen from table [l] the SU{2) correlation length Csu(2) 
and the Z(2) correlation length £z(2) show perfect agreement within errorbars. This indicates 
that the Z{2) degrees of freedom alone carry the full SU{2) MG. 

For comparison we also computed the internal energy E ( |5l| ) which corresponds to the 
correlation function of spins separated by one lattice spacing. At such a small distance and for 
large enough values of f3 the spin system is well ordered. According to a vortex condensation 
theory this means that thin vortices are suppressed and the spin configurations are not Z(2)-like 
at short distances. Thus, we do not expect agreement between the full SU(2) internal energy 
Esu(2) an d the internal energy Ez(2) extracted from Z(2) degrees of freedom. Table |2| shows 
numerical results which suggest that this is indeed the case: There is a large discrepancy between 
E§u(2) an d E z r 2 \ which indicates that the Z{2) degrees of freedom cannot account for the short 
distance behaviour of the SU (2) model. 
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L 


2(3 


%SU(2) 


£z(2) 


He 


128 


2.2 


14.01 (0.16) 


14.36 (0.24) 


14.02 (0.03) 




2.3 


18.95 (0.13) 


18.66 (0.19) 


18.91 (0.05) 




2.4 


25.39 (0.13) 


25.10 (0.18) 


25.15 (0.07) 




2.5 


33.14 (0.12) 


32.95 (0.17) 


33.19 (0.07) 




2.6 


41.77 (0.12) 


41.32 (0.19) 


41.38 (0.09) 




2.7 


49.44 (0.13) 


49.61 (0.19) 


49.60 (0.11) 


256 


2.4 


26.02 (0.37) 


25.23 (0.60) 


25.50 (0.20) 




2.5 


34.51 (0.28) 


34.63 (0.43) 


34.97 (0.16) 




2.6 


46.44 (0.25) 


46.86 (0.38) 


46.66 (0.17) 




2.7 


62.06 (0.25) 


62.28 (0.37) 


61.90 (0.23) 




2.8 


78.20 (0.26) 


78.55 (0.37) 


78.48 (0.27) 



Table 1: Estimates for the correlation length £su(2) extracted from SU(2) degrees of freedom 
and for £z(2) determined with the Z(2) degrees of freedom. The data in the right column denoted 
by £e show results for £su(2) taken from [19|. 



L 


2(5 


Esu{2) 


E 


Z(2) 


Ee 


128 


2.2 


1.24528 


(0.00027) 


0.92359 


(0.00062) 


1.24509 (0.00001) 




2.3 


1.28420 


(0.00018) 


0.95440 


(0.00044) 


1.28394 (0.00001) 




2.4 


1.31898 


(0.00016) 


0.98440 


(0.00043) 


1.31886 (0.00001) 




2.5 


1.35023 


(0.00011) 


1.01017 


(0.00042) 


1.35029 (0.00001) 




2.6 


1.37882 


(0.00009) 


1.03454 


(0.00038) 


1.37874 (0.00001) 




2.7 


1.40463 


(0.00007) 


1.05722 


(0.00042) 


1.40460 (0.00001) 


256 


2.4 


1.31541 


(0.00358) 


0.98089 


(0.00237) 






2.5 


1.34857 


(0.00171) 


1.00920 


(0.00128) 


1.35021 (0.00001) 




2.6 


1.37834 


(0.00015) 


1.03407 


(0.00034) 


1.37860 (0.00001) 




2.7 


1.40443 


(0.00006) 


1.05648 


(0.00032) 


1.40443 (0.00001) 




2.8 


1.42805 


(0.00005) 


1.07736 


(0.00031) 


1.42806 (0.00001) 



Table 2: Estimates for the internal energy Egum extracted from SU{2) degrees of freedom and 
for Ez{2) determined with the Z(2) degrees of freedom. The data in the right column denoted 
by Ee show results for Egjj^ 2 ) taken from p9fl . 
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Finally let us compare our results obtained in the 2D SU{2) principal chiral model with 
similar results obtained in AD gauge theories. In gauge theories, it was shown by numerical 
simulations that the asymptotic ST is nearly reproduced by U{\) and Z(2) degrees of freedom 
respectively. While U(l) dominance [|^] is in favour of the dual superconductor picture of 
confinement, Z(2) dominance [11] indicates that a vortex condensation mechanism is responsible 
for confinement of static quarks. However, in both cases the gauge needs to be partially fixed 
to extract the relevant degrees of freedom, and this procedure is apparently ambiguous. It 
was reported for example that the Abelian U(l) dominance is convincingly observed only in 
one particular gauge, the so-called maximal Abelian gauge. We want to emphasize that in the 
2D SU{2) principal chiral model considered in this paper, there is no gauge freedom and thus, 
there is no need to fix the gauge. The determination of the Z{2) degrees of freedom ([52]) is 
unambiguous and the same is true for the results presented in tables [l] and ^. 



4.2 Positive link model 

The numerical results of the previous section suggest that the Z{2) degrees of freedom alone 
reproduce the long distance behaviour of the SU(2) correlation function while this is not the 
case at short distances. This observation is in favour of a vortex condensation theory according 
to which the SU(2) spin field configurations behave Z(2)-like only at large distances. In this 
section we investigate the dependence of this mechanism on the short distance structure of the 
SU{2) model. We restrict the trace of the link variables to positive values and call this model 
the positive link model (PLM). According to the conventional scenario the continuum limit of 
the SU(2) model has to be taken at (3 — > oo. Thus, the PLM and the standard SU{2) model 
have the same continuum limit. Moreover, a possible mechanism which is responsible for the 
existence of a non-zero MG should be the same in both models. 

The positive link model is defined as an analogue of the positive plaquette model in gauge 
theories Q. It restricts the trace of link variables to positive values and thus suppresses thin 
vortices of the order of one lattice spacing. Such thin vortices can be seen as lattice artifacts. 
They should not influence the above discussed vortex condensation mechanism which is based 
on condensation of thick vortices having a linear extension of many lattice spacings. In this 
sense, the PLM is closer to the continuum than the standard SU(2) model. The action of the 
PLM can be written as 

Splm = ^Tr(U x l4 tn ) - A^[l- S ign(Tr(^c/t n ))], (53) 

x,n x,n 

where for a complete suppression of negative links we have to choose A = oo. To simulate the 
partition function Zplm = / Dfi(U x ) ex.p(SpLM) of the PLM we use a heatbath algorithm. In 
the update the change of a SU (2) spin variable is rejected, if the trace of one or more of the 
resulting four links is negative. Fortunately, the rejection rate decreases with increasing values of 
f3. Simulations were run on a lattice with size L = 128 and at four different j3- values. In all cases 
we started with the trivial configuration and applied 2 • 10 4 warm up sweeps. We then generated 
5 • 10 6 configurations and measured the quantities of interest in every fifth configuration. It 
should be emphasized that in this paper we are not interested in the question of scaling in the 
PLM. Our aim is to find out whether the Z{2) degrees of freedom carry the full SU(2) MG in 
the PLM as well. We thus calculated the correlation lengths Csu(2) (P]) an d £z(2) as introduced 
in the previous section. The numerical results are shown in table p[ It is clearly seen that at 
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Figure 3: Distribution of link traces in the positive link model. 



L 2/3 


£SU(2) 




E SU(2) 


E Z(2) 


128 1.7 


20.24 (0.19) 


19.17 (0.40) 


1.27350 (0.00001) 


0.93670 (0.00032) 


1.8 


23.53 (0.21) 


23.92 (0.46) 


1.29278 (0.00001) 


0.95381 (0.00043) 


1.9 


26.93 (0.27) 


27.24 (0.55) 


1.31180 (0.00001) 


0.97044 (0.00058) 


2.0 


31.63 (0.27) 


31.63 (0.94) 


1.33056 (0.00001) 


0.98686 (0.00079) 



Table 3: Estimates for the correlation length £ and the internal energy E extracted from both 
SU(2) and Z{2) degrees of freedom in the positive link model. 

least for large enough values of (3 the two quantities agree within errorbars. This shows that 
in the PLM the Z{2) degrees of freedom reproduce the long distance behaviour of the SU (2) 
correlation function as well. We interpret this result as an indication that also at weak coupling 
condensation of thick vortices is the mechanism which leads to the existence of a non-zero MG. 
For comparison we computed the internal energy Egi/M) (H)- For reasons given in the previous 
section we do not expect that the SU{2) internal energy Esu(2) an d the internal energy E^ti) 
extracted from Z{2) degrees of freedom agree. Table || shows that this is indeed the case. 

5 Summary 

In this article we discussed a vortex condensation mechanism to explain a non-zero MG in the 2D 
SU(2) principal chiral model. Following the original idea of Mack and Petkova for 4D non- 
Abelian gauge theories we formulate a sufficient condition for the MG to be non-vanishing. This 
condition is expressed in terms of the behaviour of the vortex free energy. However, the vortex 
condensation mechanism presented here does not specify the definition of a vortex. The only 
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important issue is that over some characteristic length scale spin field configurations which are 
crucial for the existence of a non-zero MG behave Z(2)-like. With this in mind we separate Z{2) 
and SO(3) degrees of freedom and construct an effective model for the Z(2) degrees of freedom 
in a SO(3) background. Assuming that Z{2) degrees of freedom alone carry the whole MG 
we arrive at an effective Ising-like expression for the SU(2) correlation function. We speculate 
that due to the non-trivial SO (3) background the effective Ising-like coupling might be smaller 
than the critical Ising coupling for all values of the bare coupling constant. In order to confirm 
the above assumption we performed numerical simulations and compared the full SU(2) MG 
with the MG extracted from Z(2) degrees of freedom. We find that Z(2) degrees of freedom 
reproduce the full SU(2) MG with perfect accuracy. This is observed not only in the standard 
model but also in the positive link model which due to the complete suppression of links with 
negative trace is closer to the continuum limit. We summarize that our numerical results as 
well as the effective model for Z{2) degrees of freedom are in favour of a vortex condensation 
mechanism. 

Finally, we want to stress the following: All arguments presented in this paper being in 
favour of a vortex theory did not use information on the phase structure of the considered 
model. Thus, they hold independently of the scenario actually being realized. In particular, 
they do not depend on whether the conventional scenario (with a non-zero MG at arbitrarily 
large (3 and asymptotic freedom) or the scenario advocated in ref. |23j (with a phase transition 
to a massless phase at a finite value of (3) is realized. Moreover, if the conjecture that the Z(2) 
MG coincides with the full SU{2) MG is correct, then our effective model might help to clarify 
this important question. 
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